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Taeyoung Lee*t, Melvin Leok*, and N. Harris McClamrochT 



o : 
o . 

(N ■ 
<D ■ 

m ; 
o\ ■ 



U 
O 



> 

in 

(N 

On 
O 
l> 
O 



Abstract — A time optimal attitude control problem is studied 
for the dynamics of a rigid body. The objective is to minimize 
the time to rotate the rigid body to a desired attitude and 
angular velocity while subject to constraints on the control 
input. Necessary conditions for optimality are developed di- 
rectly on the special orthogonal group using rotation matrices. 
They completely avoid singularities associated with local pa- 
rameterizations such as Euler angles, and they are expressed 
as compact vector equations. In addition, a discrete control 
method based on a geometric numerical integrator, referred to 
as a Lie group variational integrator, is proposed to compute 
the optimal control input. The computational approach is 
geometrically exact and numerically efficient. The proposed 
method is demonstrated by a large-angle maneuver for an 
elliptic cylinder rigid body. 

I. Introduction 

The time optimal control of spacecraft has received consis- 
tent interest as rapid attitude maneuvers are critical to various 
space missions such as military observation and satellite 
communication. The objective is to reorient the attitude of 
the spacecraft in a minimal maneuver time with constrained 
control moments. To accomplish many space missions, large- 
angle attitude maneuvering capabilities are required. 

Time optimal attitude maneuvers have been extensively 
studied in the literature [1]. The time optimal solution is 
found for a single degree of freedom system, where the atti- 
tude maneuver is constrained to an eigen-axis rotation, in [2]. 
It is known that the eigen-axis rotation is not generally time 
optimal [3], [4]. The attitude dynamics is often simplified 
in an optimality analysis, e.g., by assuming an inertially 
symmetric rigid body model [3], [4], [5], linearization [6] 
and constant magnitude angular velocity [5]. 

The attitude is defined as the orientation of a body-fixed 
frame with respect to a reference frame, and it is represented 
by a rotation matrix that lies on the special orthogonal group, 
SO(3). However, most existing optimal control scheme for 
the dynamics of a rigid body uses coordinate representations 
such as Euler angles and quaternions. The minimal attitude 
representations like Euler angles and Rodrigues parame- 
ters have singularities, so they are not desirable for large- 
angle maneuvers. The non-minimal attitude representations 
like quaternions have associated problems. Besides the unit 
norm constraint, the quaternion representation double covers 
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SO(3). So, it has an inevitable ambiguity in expressing the 
attitude. 

The objective of this paper is to solve the time optimal 
attitude control problem directly on SO (3) using rotation 
matrices without need of any attitude parameterization. Using 
a specific property of the special orthogonal group, namely 
that the Lie algebra so(3) is isomorphic to K 3 , necessary con- 
ditions for optimality are developed and represented as vector 
equations on R 3 . They avoid singularities associated with 
Euler angles completely, and the resulting expressions for 
the optimality necessary conditions are more compact than 
expressions obtained by using quaternions. Consequently, the 
attitude dynamics need not be simplified to make the optimal 
control problem tractable. 

The remaining part of this paper is focused on devel- 
oping a computational approach to solve this optimal con- 
trol problem. The dynamics of a rigid body has certain 
geometric features; in addition to the configuration space 
being a Lie group, the dynamics are characterized by sym- 
plectic, momentum and energy preserving properties. The 
most common numerical integration methods, including the 
widely used (non-symplectic) explicit Runge-Kutta schemes, 
preserve neither the Lie group structure nor these geometric 
properties. 

Lie group variational integrators are geometric numerical 
integrators that preserve these geometric features of the 
rigid body dynamics [7]. Based on this structure-preserving 
numerical integrator, computational approaches have been 
proposed to solve various optimal control problems for the 
dynamics of rigid bodies [8], [9], [10]. In this paper, the time 
optimal attitude control problem is discretized at the level 
of the initial problem formulation, and discrete necessary 
conditions for optimality are developed using the Lie group 
variational integrator. This provides geometrically exact but 
computationally efficient tools. 

In summary, the optimization scheme for time optimal 
attitude maneuvers that we present in this paper has the 
following important features: (i) necessary conditions for 
optimality are developed directly on SO(3), and (ii) a compu- 
tational approach is adopted by using a Lie group variational 
integrator for overall numerical accuracy and efficiency. 

This paper is organized as follows. The time optimal 
attitude control problem is formulated, and continuous-time 
necessary conditions for optimality are developed in Section 
II, and in a parallel fashion, a discrete-time optimal control 
method is presented in Section III, followed by numerical 
examples in Section IV. 



II. Time Optimal Attitude Control 

A. Equations of Motion 

We consider the attitude dynamics of a rigid body. The 
configuration space is the special orthogonal group SO (3), 

SO(3) = {R £ R 3x3 | R T R = I 3x3 , deiR=l}, 

where the rotation matrix R 6 SO (3) represents the linear 
transformation from the body-fixed frame to the inertial 
frame. 

The continuous equations of motion for the attitude dy- 
namics of a rigid body are given by 



jfi + fi x Jfi = u, 
R = R&, 



(1) 

(2) 



where the matrix J £ R 3x3 is the moment of inertia matrix, 
the vector fi £ R 3 is the angular velocity expressed in the 
body-fixed frame, and the external control moment is denoted 
by u £ R 3 . The hat map ~ : R 3 i— > so(3) is an isomorphism 
from R 3 to skew-symmetric matrices so(3), and is defined 
by the condition xy = x X y for all x, y £ R 3 . The inverse 
map is denoted by the vee map (-) v : so (3) i— » R 3 . 

B. Time Optimal Attitude Control Problem 

The objective of the time optimal attitude control problem 
is to transfer the given initial attitude and the angular velocity 
(R 0) Cl ) of the rigid body to the desired values (Rf,£lf) 
within a minimal maneuver time t f with constrained control 
moment ||u|| < u for a given control limit u £ R. 

Forgiven: (R a , fi D ), (R f , fi/), u 
min <J=l Idt 



such that R(t } ) = R f , fi(t/) = fi/, 
subject to \\u(t)\\ <u Vt e [0,t/] and JTJ, ©. 

C. Necessary Conditions for Optimality 

We solve this optimal control problem using variational 
principles applied on SO(3). Expressions for variations of a 
rotation matrix, and transversality conditions are presented, 
and necessary conditions for optimality are developed. 

Expressions for variations: We represent a variation of 
a rotation matrix using the exponential map, cxp : so (3) i— > 
SO (3) 



R e = i?exp ef), 



(3) 



where e £ (— c, c) for c > 0, and f) £ so(3) for r\ £ R 3 . 
Since the exponential map is a local diffeomorphism, this 
expression is well-defined for some constant c for given f). 
The infinitesimal variation of the rotation matrix is given by 

d 



SR 



de 



Rexpefj = Rfj. 



(4) 



-=o 



The infinitesimal variation of R T R is obtained from (O and 
© as 

5(R T R) = 5R T R + R T 5R, 



= -r]R T R + R t {Rt) + Rrj), 
= f) + flrj — f)f2, 

= (77 + Q x rtf. (5) 

The variational expressions given by (O and (0) are the key 
ingredients to developing necessary conditions for optimality 
for an arbitrary optimal attitude maneuver. 

Transversality conditions: The differentials in the termi- 
nal attitude and the terminal angular velocity are composed 
of the variation for a fixed time and a term due to the terminal 
time variation. Since the terminal boundary conditions are 
fixed, we have the transversality conditions as 

SR(t f ) + R(t f )dt f = R(t f )fj(t f ) + R(t f )dt f = 0, (6) 
5ft(t f ) + (l(t f )dtf = 0. (7) 

Necessary conditions for optimality: Define the aug- 
mented performance index as 

rtf 



J a = 1 + X ■ (u — fl X Jfl — JQ) 

Jo 

+ \ R ■ (Cl- R T R) y dt, 

where X n , X R £ R 3 are Lagrange multipliers. 

Using the infinitesimal variation of the augmented 
performance index is given by 

sj a = x n ■ {Su -ifOxjn-flx jsn - jsn) 

Jo 

+ X R ■ • {SQ - fi - Q x 77) dt 

+ {1 + a° • (u - n x jn - jo) + x B ' ■ (n - i? T i?) v } dtf. 

tf 

Using integration by parts, we obtain 

sj a = x n ■ (5u - sn x jo - fi x jsn) + x n ■ jsn 

Jo 

+ x R -(sn-nx^ + x 11 -r]dt 

- {x n ■ jsn + x R ■ t}} tf 


+ {l + X n ■ (u - n x Jfi - Jfi) + X R ■ (Cl - R T R) V } dtf. 

tf 

Since the initial attitude and the initial angular velocity 
are fixed, we have 77(0) = 0, Jfi(0) = 0. Substituting 
and rearranging, the infinitesimal variation of the augmented 
performance index is given by 

SJ a = I 1 SSI- {-Jfi x A° - J(X n x fi) + JX n + X R } 



+ 77 • jfi x X R + A fl | + 5u ■ X n dt 

+ {l + X n ■ (u - fi x Jfi) + X R ■ fi} dtf. 

tf 

We choose multiplier equations and boundary conditions 
such that the expressions in all braces in the above equations 
are identically zero. Then, we have 



SJa 



Su ■ X n dt. 



The optimal control input u must satisfy 

A n • Su > 0, 



(8) 



for all admissible Su in t 6 [0, t/]. If X n = for a finite 
time period, the control input is not determined by {8). Such 
solutions are referred to as singular arcs. Later, it is shown 
that there is no singular arc in this optimal control problem. 

In summary, the necessary conditions for optimality are 
given by 
• Multiplier equations 



JX n + J(fi x A n ) - jn x X Q + X R 
X R + x X R = 0, 

• Optimality condition 

u = -u(X n /\\X n \\), 

• Boundary and transversality conditions 

OR(0),fi(0)) - (iZo.fio), 

(ii(i / ),n(t / )) = ( J R / ,fi / ), 

{1 + A" • (u-Cl x JQ) + A fi • 0} 



0. 



= 0, 



(9) 
(10) 

(ID 

(12) 
(13) 
(14) 



Assuming that the rigid body is inertially symmetric, J = 
^3x3, the multiplier equation (0 is reduced to A° + X R = 0. 

These necessary conditions for optimality are valid for atti- 
tude maneuvers of arbitrary magnitude as they are developed 
by using the rotation matrix representation on SO (3). Since 
the variation of the rotation matrix is expressed in terms 
of the Lie algebra so (3) isomorphic to R 3 , the multiplier 
equations are written as compact vector equations on M 3 . The 
presented necessary conditions for optimality have neither 
the singularities inherent to Euler angles nor the ambiguities 
and redundancy associated with quaternions. 

D. Singular arc 

In this subsection, we show that singular arcs do not 
exist along a solution of this time optimal control problem. 
Suppose that there exist a singular interval, i.e. A (t) = for 
a finite time period in [0, tt]. Then, the minimum principle 
given by © does not lead to a well-defined condition for the 
optimal control input. Instead, the control input is determined 
by the requirement that the time derivative of A is equal to 
zero. 

Let the 2g-th time derivative of X n be the lowest order 
derivative in which the control input u appears explicitly 
with a coefficient that is not identically zero on the singular 
interval. Then, the integer q is called the order of the singular 
arc [11]. Here, due to the special linear structure of this 
multiplier equation, the singular arc has infinite order. If the 
condition X n = A = is satisfied at a single point along the 
trajectory, X R = X R = 0, and these are satisfied identically 
throughout the trajectory independent of the control input. In 
this case, it is clear that the boundary condition ([141 cannot 
be satisfied. Thus, there is no singular arc in an optimal 
solution. 



III. Discrete-time Time Optimal Attitude 
Control 

In this section, we present a computational approach, 
referred to as discrete optimal control of discrete Lagrangian 
systems [12], to solve the time optimal attitude control prob- 
lem numerically. In this approach, the dynamics of the rigid 
body is discretized using the discrete Hamilton's principle, 
in order to obtain a Lie group variational integrator [7]. The 
corresponding discrete equations of motion are imposed as 
dynamic constraints to be satisfied by using Lagrange multi- 
pliers, and necessary conditions for optimality, expressed as 
discrete equations on multipliers, are obtained. 

This method yields substantial computational advantages 
in finding an optimal control solution. The discrete dynamics 
are more faithful to the continuous equations of motion, and 
consequently more accurate solutions to the optimal control 
problem are obtained. It has been shown that the discrete 
dynamics is more reliable even for controlled system as it 
computes the energy dissipation rate of controlled systems 
more accurately [13]. In particular, the discrete flow of the 
Lie group variational integrator remains on SO (3). 

Optimal solutions, computed using an indirect approach, 
are usually sensitive to small variations of the multipliers. 
This causes difficulties, such as numerical ill-conditioning, 
when solving the necessary conditions for optimality ex- 
pressed as a two-point boundary value problem. Sensitivity 
derivatives, computed using the discrete necessary condi- 
tions, are not corrupted by numerical dissipation caused by 
conventional numerical integration schemes. Thus, the pro- 
posed computational approach is more numerically robust, 
and the necessary conditions can be solved in a computa- 
tionally efficient manner. 

A. Lie Group Variational Integrator 

Since the dynamics of a rigid body has the structure of 
a Lagrangian or Hamiltonian system, they are symplectic, 
momentum and energy preserving. These geometric features 
determine the qualitative behavior of the rigid body dynam- 
ics, and they can serve as a basis for theoretical study of 
rigid body dynamics. 

In contrast, the most common numerical integration meth- 
ods, including the widely used (non-symplectic) explicit 
Runge-Kutta schemes, preserve neither the Lie group struc- 
ture nor these geometric properties. Additionally, if we 
integrate <(2J using a typical Runge-Kutta scheme, the quan- 
tity R T R inevitably drifts from the identity matrix as the 
simulation time increases. 

In [7], Lie group variational integrators are constructed by 
explicitly adapting Lie group methods [14] to the discrete 
variational principle [13]. They have the desirable property 
that they are symplectic and momentum preserving, and they 
exhibit good energy behavior for an exponentially long time 
period. They also preserve the Lie group structure without 
the use of local charts, reprojection, or constraints. These ge- 
ometrically exact numerical integration methods yield highly 
efficient and accurate computational algorithms for rigid 
body dynamics, and avoid singularities and ambiguities. 



Using the results presented in [7], a Lie group variational 
integrator on SO (3) for equations (|T), © is given by 

hJSl k = F k J d - J d F£, (15) 
Rk+i = R k F k , (16) 
,m k+1 = F^JQk + huk+x, (17) 

where the subscript k denotes the £>fh step for a fixed 
integration step size h € K. The matrix Jd € M 3x3 is 
a nonstandard moment of inertia matrix defined by = 
itr[J]7 3x3 - J 6 M 3x3 . The matrix F h G SO(3) denotes 
the relative attitude between adjacent integration steps. 

For given (R k ,x k ) and control input, $15[ is solved to 
find F k . Then (R k+ i, Qk+i) are obtained by (fTol i and ( fTTI i. 
This yields a map (R k ,fi k ) >— > (-Rfc+i, Ofc+i), and this 
process is repeated. The only implicit part is (fl3T l. where the 
actual computation of F k is done in the Lie algebra so (3) 
of dimension 3. 

One of the distinct features of the Lie group variational 
integrator is that it preserves both the symplectic property 
and the Lie group structure of the rigid body dynamics. 
As such, it exhibits substantially improved computational 
accuracy and efficiency compared with other geometric inte- 
grators that preserve only one of these properties such as non- 
symplectic Lie group methods [15]. The symplectic property 
is important even in the case of controlled dynamics, since 
the dissipation rate of the total energy is typically computed 
inaccurately by non-symplectic integrators [13]. 

B. Discrete-time Time Optimal Attitude Control Problem 

The objective is to transfer the rigid body in a prescribed 
way within a minimal discrete maneuver time N with 
constrained control input. 

Forgiven: (R ,Q, ),(R f ,Q.f),u 

mini J = V l), 

such that Rn = Rf, H.^ = {If, 
subject to |K + i|| < u Vfc G [0,iV- 1] and (|T5 ) -(fT7 ]) , 

C. Discrete-Time Necessary Conditions for Optimality 

Expressions for variations: Similar to |@), the variation 
of rotation matrices Rk and Fk are expressed as 

5R k = R k fi k , SF k = F k i k (18) 

for r]k,£,k S K 3 . Using this and (1161 1. the variation of 
RjRk+i is given by 

S(R k R k+ i) — SR^Rk+i + R^SRk+i, 
= -vFk + F k rj k+1 , 

= F fe (-if % + %+1 r, d9) 

where the property F T x = F T xF for any x G K 3 and 
F G SO (3) is used in the last step. 

Now we develop an expression for a constrained variation 
corresponding (I151 l. Taking a variation of (Q3), we obtain 

hjm~ k = F k i k jd + JdikF^. 



Using the property, xA + A T x = ({tr[A] ^3x3 — A} xf~ for 
all x G K 3 A G M 3x3 , the above equation can be written as 

hJSd k = FkTkFkJd + JdFfFkTk, 

= (MF k Jd] 7 3 x3 - F k Jd}F k ^ k y. 
Thus, the vector £/. is expressed in terms of 5Vt k 

a = B k JSn k , (20) 

where B k = hF? {tr[F fe J d ] 7 3 x3 - F^y 1 G K 3x3 . This 
shows the relationship between 5il k and 5F k . 

Transversality conditions: Similar to (O, we choose the 
transversality conditions for the angular velocity as 

5Q N + (fijy - 0;v-i)<W = 0. (21) 

The variation of the terminal attitude due to the terminal time 
change is expressed as 

Rn I 2^-n-i(Rn — Rn-i) + -^R-n(Rn — Rn-i) J- 5N 

= ^Rn{F N -i~FT_ 1 }6N. 

This expression is chosen such that it respects the skew- 
symmetry of a Lie algebra so(3) element. Using this, the 
transversality conditions for the attitude are given by 

R n Vn + ^Rn{Fn-i-FJ;_ 1 }SN = 0, (22) 

Necessary conditions for optimality: Define the aug- 
mented performance index as 

N-l 

Ja=^l + X k - {-J^k+1 + F^Jflk + hu k+1 } 
k=Q 

+ A? • ~ ((F k F?y (RT Rk+1 R T k+1 R k y) . 

Here we assume that the time step size h is small so that the 
relative attitude rotation between adjacent integration steps 
is less than i.e. ||(logmFfc) v || < Then, F k is equal to 
R k R k +i if and only if their skew parts are identical, which 
can be easily shown using Rodrigues' formula. Equation ( fTBT l 
is considered implicitly using a constrained variation. 

Using (T% , the infinitesimal variation of the augmented 
performance index is given by 

JV-l 

sj a = J2 x k- i hSu k+i - Jsn k+1 + 5F^m k + F^jsrik} 

k=0 

+ X^l{F k ^ k + F k T Vk ^r lk+1 r 

+ (a+F k T Vk - Vk+1 rF k T Y 

+ {1 + A^_ : • {-JQ N + F^_ t Jfijv-i + hu N }}SN 
+ X%-i-l(F N -i-F%_ 1 ) v 5N 

~ ^n-i ' ^Rn-iRn — RJjRn-i) w SN. 

Several algebraic manipulation steps are required here; (i) 
using the property xA+A T x = {{ti[A] I 3x3 — A} xf for all 



i£l 3 and A £ R 3x3 , the expression in the second braces 
is written as a vector form, (ii) equation ( f20b is substituted 
to express in terms of 6&l k , and (iii) using the fact that 
r/o = 0, (50o = 0, the summation indices for the variables at 
the k + 1-th step are rewritten, which is the discrete analog 
of integration by parts. Then, we obtain 

N-l 

5J a = X k ■ h5u k+1 

k=Q 

N-l 

+ £ «j fc . { - + J(F k - BlFfm k )\% 

k=\ 

+ l -JB T k (^F k ]I-F k )\^) 

N-l , 

+ E%-{2 (tr[Ffc - i]/ " Ffc - i)A ^ Li 

k=l 

-l Fk (ti{F k }I-F k )\x} 

- X N _ t ■ JSQ N - A£_! ■ ^HF N ^} I - F^JriN 
+ {1 + A^_! • {- JQ N + F^_ 1 jn N - 1 + hu N }}6N 
+ \ N - 1 -\(F N - 1 -F?_ i ySN 

~ ^n-i ' ^{Rn-iRn — R n Rn-x) v 6N. (23) 

Substituting the transversality conditions (f2TT > and ( l22l . all 
of the expressions in the last four lines of the above equation 
are reduced to 

{l + A^_! • {- Jftjv-i + F§_ x JOjv-i + hu N ) 

+ *n-i ■ \ ((Fn-i) 2 - (F^ff }SN. (24) 

We choose discrete multiplier equations such that the 
expressions in the first two braces in (1231 are identically zero, 
and we choose boundary condition such that the expression 
given by (l24l is equal to zero. Then, we have 

N-l 

SJa = Yl A fc ' h6u k+l- 

k=0 

The optimal control input u k +\ must satisfy 

• Su k+1 > 0, 

for all admissible 5u k +i and k G {0, • • ■ , N — 1}. Here, we 
do not show that there is no singular arc in the discrete-time 
optimal control problem. We assume that the result presented 
in Section III-DI for the continuous-time case also applies to 
the discrete-time case. In summary, the discrete necessary 
conditions for optimality are given by 

• Multiplier equations 

- JAti + J(F k - BlFfm k )\% 

1 (25) 
+ -JBl(tiF k ]I-F k )\f = 0, 

(tr[F fc _!] / - fWijAf.! - F fc (tr[F fc ] / - F k )\« = 0. (26) 



• Optimality condition 

u k+1 = -u(\2/\\\%\\) (27) 

• Boundary and transversality conditions 

(iio.no) = (RoM (28) 

(R N ,n N ) = (R f> n f ), (29) 

1 + X N-i ■ {-J^N-i + i^_i J^jv-i + hu N ) 
+ *N-i-\{{F N _ 1 f-{FT_ l ?) y =Q. (30) 

In the above equations, the only implicit part is ( TT3T ). For a 
given initial condition {(iio, Oo)j (Aq , Aq )}, we solve ( fT~5T > 
to obtain Fq, and we find the control input u\ by d27l i. 
Then, (i?i,Qi) are obtained by ( TT~6b and (fTTI i. Using f2i, 
we solve (fT~5b to obtain Fi. Finally, (A*, A^ 2 ) are obtained 
by 65) and (g5). This yields a map {(i? , n ), (X§, Aq )} ^ 
{(i?i,J7i), (A^, A^)}, and this process is repeated. 

The discrete necessary conditions for optimality are given 
by a two-point boundary value problem. This is to find the 
optimal discrete flow, multiplier, control input, and terminal 
maneuver time to satisfy the equations of motion (fT5ll-(fT7]l. 
multiplier equations (l25T l, (|26| |. optimality condition ( |27] |. and 
boundary conditions (T28b — (f30b simultaneously. 

We use a neighboring extremal computational 
method [16]. A nominal solution satisfying all of the 
necessary conditions except the boundary conditions is 
chosen. The unspecified initial multiplier is updated so as 
to satisfy the specified terminal boundary conditions in 
the limit. This is also referred to as a shooting method. 
The main advantage of the neighboring extremal method 
is that the number of iteration variables is small. In 
other approaches, the initial guess of control input history 
or multiplier variables are iterated, so the number of 
optimization parameters are proportional to the number of 
discrete time steps. 

A difficulty is that the extremal solutions are sensitive 
to small changes in the unspecified initial multiplier values. 
The nonlinearities also make it hard to construct an accurate 
estimate of sensitivity, and it may result in numerical ill- 
conditioning. By adopting a geometric numerical integrator, 
sensitivity derivatives along the discrete necessary conditions 
do not have numerical dissipation introduced by conventional 
numerical integration schemes. Thus, they are numerically 
more robust, and the necessary conditions can be solved 
computationally efficiently. 

IV. Numerical Example 

We choose an elliptic cylinder for a rigid body model 
with semi-major axis 0.8 m, semi-minor axis 0.2 m, height 
0.6m, mass l,kg. The moment of inertia matrix is J = 
diag[0.04, 0.19, 0.17] kgm 2 , and the maximum control in- 
puts is chosen as u — 0.1 Nm. 

The desired attitude maneuver is a rest-to-rest large angle 
rotation given by 

(flo.no) = (isxa.O) 
(Rf,Sl f ) = (expfo.O), 



(a) Attitude maneuver (a) Attitude maneuver 




(d) Lagrange multiplier A n (e) Lagrange multiplier \ R (d) Lagrange multiplier A n (e) Lagrange multiplier \ R 

Fig. 1. Time optimal attitude maneuver, 8 = 120° Fig. 2. Time optimal attitude maneuver, 8 = 180° 



where v = -75 [lj 1, 1] G K 3 , and is varied as 120° and 
180°. 

When deriving the discrete necessary conditions for opti- 
mality, we assume that the number of discrete steps N varies. 
For computational purpose, it is not desirable to search the 
optimal value of N since the terminal attitude, angular ve- 
locity and multiplier change in a discrete manner for varying 
integer N. Thus, it is not guaranteed that the boundary 
condition is satisfied to a desired numerical accuracy. 

In the numerical computation, we fix the number of steps 
by an educated guess, N = 1000 in this particular numerical 
example, and we vary the timestep h. In essence, we find the 
seven parameters, initial multiplier (A^, Aq) and the time 
step h, satisfying the seven-dimensional terminal boundary 
conditions (I28ll-(l30ll under the discrete equations of motion, 
the multiplier equation, and the optimality condition. 

We solve this two-point boundary value problem, inter- 
preted as a nonlinear equation by the shooting method, 
using a general nonlinear equation solver, namely the Matlab 
f solve function. The multipliers are initialized randomly, 
and the timestep is initialized as h = 0.002 seconds. 
The optimal solutions are found in 94 and 211 seconds, 
respectively, on Intel Pentinum M 1.73 GHz processor, and 
the boundary condition errors are less than 10 -15 . 

The optimized attitude maneuver, angular velocity, 



multiplier, and control input history are presented 
in Figures Q] and [2] (Simple animations which show 
these maneuvers of the rigid body are available at 
|http : / /www . umich . edu/~tylee| ) The optimized 
maneuver times are 3.3855 and 3.8184 seconds, respectively, 
and there is no singular arc along the optimized solutions. 

V. Conclusions 

A time optimal attitude control problem to rotate a rigid 
body within a minimal time with constrained control input 
is studied. Necessary conditions for optimality are developed 
on SO (3) using rotation matrices without need of attitude 
parameterizations such as Euler angles and quaternions. This 
provides a globally applicable and compact form of nec- 
essary conditions for optimality. For overall computational 
accuracy and efficiency, a discrete optimal control method is 
proposed using a Lie group variational integrator. 

In this paper, the two-norm of the control moment is 
constrained, and consequently, there is no singular arc in 
the optimal solution. The proposed necessary conditions for 
optimality can be directly applied, without modification, to 
the case where the absolute value of each component of the 
control moment is bounded. In this case, the expressions 
for optimal singular control can be developed, for example, 
by following the approach given in [4], using the compact 
multiplier equations presented here. 
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